subroutine init_sink_data(lunit)
  use data_param,only:geom_exist
  use data_param,only:sink_exist
  use label_data_param,only:n_obj_type
  use geom_data_param,only:n_geom_obj
  use sink_data_param,only:n_sink_type
  use sink_data_param,only:sink_type
  use sink_data_param,only:sink_map
  implicit none

  integer,intent(in) :: lunit
  integer :: i,j,ierr

  n_sink_type=0
  nullify(sink_type)
  nullify(sink_map)

  if (.not.geom_exist) return
!  if (.not.sink_exist) return
  if (n_geom_obj.eq.0) return

  allocate(sink_map(n_geom_obj,n_obj_type),stat=ierr)

  if (ierr.ne.0) call abnormal_stop("init_sink_data",&
       "Memory allocation causes an error.")

  do i=1,n_geom_obj
     do j=1,n_obj_type
        sink_map(i,j)=0
     end do
  end do

  return

end subroutine init_sink_data
